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Abstract 

The braneworld description of our universe entails a large extra dimension and a fundamental 
scale of gravity that might be lower by several orders of magnitude as compared to the Planck scale. 
An interesting consequence of the braneworld scenario is in the nature of the vacuum solutions of 
the brane gravitational field equations, with properties quite distinct as compared to the standard 
black hole solutions of general relativity. One possibility of observationally discriminating between 
different types of black holes is the study of the emission properties of the accretion disks. In the 
present paper we obtain the energy flux, the emission spectrum and accretion efficiency from the 
accretion disks around several classes of static and rotating brane world black holes, and we compare 
them to the general relativistic case. Particular signatures can appear in the electromagnetic 
spectrum, thus leading to the possibility of directly testing extra-dimensional physical models by 
using astrophysical observations of the emission spectra from accretion disks. 
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I. INTRODUCTION 



The idea, proposed in l|,|2j, that our four- dimensional Universe might be a three-brane, 
embedded in a five-dimensional space-time (the bulk), has attracted a considerable interest 
in the past few years. According to the brane-world scenario, the physical fields (electro- 
magnetic, Yang-Mills etc.) in our four- dimensional Universe are confined to the three brane. 
These fields are assumed to arise as fluctuations of branes in string theories. Only gravity can 
freely propagate in both the brane and bulk space-times, with the gravitational self-couplings 
not significantly modified. This model originated from the study of a single 3-brane embed- 
ded in five dimensions, with the 5D metric given by ds 2 = e~^ y ^r\ llv dx lL dx v + dy 2 , which, 
due to the appearance of the warp factor, could produce a large hierarchy between the scale 
of particle physics and gravity. Even if the fifth dimension is uncompactified, standard 
4D gravity is reproduced on the brane. Hence this model allows the presence of large, or 
even infinite non-compact extra dimensions. Our brane is identified to a domain wall in a 
5-dimensional anti-de Sitter space-time. For a review of the brane world models see [31]. 

The braneworld description of our universe entails a large extra dimension and a funda- 
mental scale of gravity that might be lower by several orders of magnitude compared to the 
Planck scale [lluj]. Due to the correction terms coming from the extra dimensions, signifi- 
cant deviations from the Einstein theory occur in brane world models at very high energies 
{4, 5]. Gravity is largely modified at the electro- weak scale of 1 TeV. The cosmological and 
astrophysical implications of the brane world theories have been extensively investigated in 
the physical literature 

Several classes of spherically symmetric solutions of the static gravitational field equations 



in the vacuum on the brane have been obtained in 

ay , [h], As a possible physical 
application of these solutions the behavior of the angular velocity v tg of the test particles 
in stable circular orbits has been considered 0, [h], . The general form of the solution, 
together with two constants of integration, uniquely determines the rotational velocity of 
the particle. In the limit of large radial distances, and for a particular set of values of 
the integration constants the angular velocity tends to a constant value. This behavior is 
typical for massive particles (hydrogen clouds) outside galaxies, and is usually explained by 
postulating the existence of the dark matter. The exact galactic metric, the dark radiation, 
the dark pressure and the lensing in the flat rotation curves region in the brane world scenario 
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has been obtained in 



111. 



For standard general relativistic spherical compact objects the exterior space-time is 
described by the Schwarzschild metric. In the five dimensional brane world models, the high 
energy corrections to the energy density, together with the Weyl stresses from bulk gravitons, 



imply t 
metric 



rat on the brane the exterior metric of a static star is no longer the Schwarzschild 



121 ] . The presence of the Weyl stresses also means that the matching conditions do 



not have a unique solution on the brane; the knowledge of the five-dimensional Weyl tensor 
is needed as a minimum condition for uniqueness. 

It is known that the Einstein field equations in five dimensions admit more general spher- 
ically symmetric black holes on the brane than four- dimensional general relativity. Hence 
an interesting consequence of the braneworld scenario is in the nature of the spherically 
symmetric vacuum solutions to the brane gravitational field equations, which could repre- 
sent black holes with properties quite distinct as compared to ordinary black holes in four 
dimensions. Such black holes are likely to have very diverse cosmological and astrophysi- 
cal signatures. Static, spherically s ym metric exterior vacuum solutions of the brane world 
models have been proposed first in [12] and in [13J]. The first of these solutions, obtained in 



12J , has the mathematical form of the Reissner- Nordstrom solution of the standard general 



relativity, in which a tidal Weyl parameter plays the role of the electric charge of the general 
relativistic solution. The solution has been obtained by imposing the null energy condition 
on the 3-brane for a bulk having non zero Weyl curvature, and it can be matched to the 
interior solution corresponding to a constant density brane world star. A second exterior 
solution, which also matches a constant density interior, has been derived in 

Two families of analytic solutions of the spherically symmetric vacuum brane world model 
equations (with g u 7^ — l/g rr ), parameterized by the ADM mass and a PPN parameter j3 
have been obtained in [14|. Non-singular black-hole solutions in the brane world model 
have been considered in [151 ] . by relaxing the condition of the zero scalar curvature but 
retaining the null energy condition. The four- dimensional Gauss and Codazzi equations for 
an arbitrary static spherically symmetric star in a Randall-Sundrum type II brane world 
have been completely solved on the brane in [lfjj. The on-brane boundary can be used to 
determine the full 5-dimensional space-time geometry. The procedure can be generalized to 
solid objects such as planets. 

A method to extend into the bulk asymptotically flat static spherically symmetric brane- 
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world metrics has been proposed in 17j. The exact integration of the field equations along 
the fifth coordinate was done by using the multipole (1/r) expansion. The results show that 
the shape of the horizon of the brane black hole solutions is very likely a flat "pancake" for 
astrophysical sources. 

The general solution to the trace of the 4-dimensional Einstein equations for static, spher- 
ically symmetric configurations has been used as a basis for finding a general class of black 
hole metrics, containing one arbitrary function g u = A(r), which vanishes at some r = > 



(the horizon radius) in 18 



Under certain reasonable restrictions, black hole metrics are 
found, with or without matter. Depending on the boundary conditions the metrics can 
be asymptotically flat, or have any other prescribed asymptotic. The exact stationary and 
axisymmetric solutions describing charged rotating black holes localized on a 3-brane in the 
Randall- Sundrum braneworld were studied in [19] . By taking the metric on the brane to be 
of the Kerr-Schild form it can be shown that the Kerr-Newman solution of ordinary general 
relativity in which the electric charge is superseded by a tidal charge satisfies a closed system 
of the effective gravitational field equations on the brane. The negative tidal charge may 
provide a mechanism for spinning up the black hole so that its rotation parameter exceeds 
its mass. For a review of the black hole properties and of the lensing in the brane world 
models see j^ . 

It is generally expected that most of the astrophysical objects grow substantially in mass 
via accretion. Recent observations suggest that around most of the active galactic nuclei 
(AGN's) or black hole candidates there exist gas clouds surrounding the central compact 
object, and an associated accretion disc, on a variety of scales from a tenth of a parsec to a 
few hundred parsecs 21|. These clouds are assumed to form a geometrically and optically 
thick torus (or warped disc), which absorbs most of the ultraviolet radiation and the soft 
X-rays. The gas exists in either the molecular or the atomic phase. The most powerful 
evidence for the existence of super massive black holes comes from the very-long baseline 
interferometry (VLBI) imaging of molecular H2O masers in the active galaxy NGC 4258 
22] . This imaging, produced by Doppler shift measurements assuming Keplerian motion 
of the masering source, has allowed a quite accurate estimation of the central mass, which 
has been found to be a 3.6 x 1O 7 M super massive dark object, within 0.13 parsecs. Hence, 
important astrophysical information can be obtained from the observation of the motion of 
the gas streams in the gravitational field of compact objects. 
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The determination of the accretion rate for an astrophysical object can give a strong 
evidence for the existence of a surface of the object. A model in which Sgr A*, the 3.7 x 
1O 6 M super massive black hole candidate at the Galactic center, may be a compact object 
with a thermally emitting surface was considered in 2s|. For very compact surfaces within 
the photon orbit, the thermal assumption is likely to be a good approximation because 
of the large number of rays that are strongly gravitationally lensed back onto the surface. 
Given the very low quiescent luminosity of Sgr A* in the near-infrared, the existence of a 
hard surface, even in the limit in which the radius approaches the horizon, places a severe 
constraint on the steady mass accretion rate onto the source, M < 1O _12 M yr~ x . This 
limit is well below the minimum accretion rate needed to power the observed submillimeter 
luminosity of Sgr A*, M > l(r lo M yr. 

Thus, from the determination of the accretion rate it follows that Sgr A* does not have a 
surface, that is, it must have an event horizon. Therefore the study of the accretion processes 
by compact objects is a powerful indicator of their physical nature. 

The first comprehensive theory of accretion disks around black holes was constructed in 



24j . This theory was extended to the general relativistic models of the mass accretion onto 



rotating black holes in [251 ] . These pioneering works developed thin steady-state accretion 
disks, where the accreting matter moves in Keplerian orbits. The hydrodynamical equilib- 
rium in the disk is maintained by an efficient cooling mechanism via radiation transport. 
The photon flux emitted by the disk surface was studied under the assumption that the disk 
emits a black body radiation. The properties of radiant energy flux over the thin accretion 
disks were further analyzed in |26| and in j^, where the effects of the photon capture by 
the hole on the spin evolution were presented as well. In these works the efficiency with 
which black holes convert rest mass into outgoing radiation in the accretion process was also 
computed. 

The emissivity properties of the accretion disks have also been investigated for exotic 



central objects recent 



rotating cases [28 



29 



v such as quark, boson or fermion stars for both rotating and non- 



301, as well as for the modified f(R) type theories of gravity [3jJ. The 



radiation power per unit area, the temperature of the disk and the spectrum of the emitted 
radiation were given, and compared with the case of a Schwarzschild black hole of an equal 
mass. 

It is the purpose of the present paper to study the matter accretion by brane world black 
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holes. By using the general formalism of accretion we analyze the accretion process for 
several black hole type solutions of the gravitational field equations on the brane, both non- 
rotating and rotating, which have been previously obtained. Particular signatures can appear 
in the electromagnetic spectrum, thus leading to the possibility of directly testing extra- 
dimensional physical models by using astrophysical observations of the emission spectra 
from accretion disks. 

The present paper is organized as follows. We review the field equations of the brane 
world models and the static, spherically symmetric solutions of the field equations as well as 
the rotating ones in Section II. The thin accretion disks onto black holes are briefly described 
in Section III. In Section IV we consider the radiation flux, spectrum and efficiency of thin 
accretion disks onto several classes of brane world black holes. We discuss and conclude our 
results in Section V. 

II. THE GRAVITATIONAL FIELD EQUATIONS IN THE BRANE WORLD 
MODELS 

In the present Section we briefly describe the basic mathematical formalism of the brane 
world models, and we present the spherically symmetric static vacuum field equations. The 
solutions of the vacuum field equations on the brane physically describe the brane world 
black holes. 

A. The gravitational field equations on the brane 

We start by considering a five dimensional (5D) spacetime (the bulk), with a single four- 
dimensional (4D) brane, on which matter is confined. The 4D brane world (^M,g^ v ) is 



located at a hypersurface [B \X A j = Oj in the 5D bulk spacetime (^M, Qab), of which 

coordinates are described by X A , A = 0, 1, ...,4. The induced 4D coordinates on the brane 
are = 0,1,2,3. In the present paper the capital Latin indices A, B, I, J, ... take 

values in the range 0, 1, 4, while the Greek indices run in the range 0, 3. 
The action of the system is given by S = Sbuik + Sbrane, where 
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2kl 



6 



and 



brane 



< 4 )M 



-(5), 



K ± + L 



brane 



d x, 



(1) 



where k\ = 87rG 5 is the 5D gravitational constant, and ^L m are the 5D scalar curvature 
and the matter Lagrangian in the bulk, Lb ra ne (#0/3, is the 4D Lagrangian, which is given 
by a generic functional of the brane metric g a p and of the matter fields ip, K is the trace of 
the extrinsic curvature on either side of the brane, and A5 and A;, (the constant brane tension) 
are the negative vacuum energy densities in the bulk and on the brane, respectively 4|. 
The Einstein field equations in the bulk are given by 



(5) G 



1. 1 



kf^T, 



u, 



4] 



^T u = -A^gu + 5(B) \-X b {5) gij + T u 



(2) 



where ^T u = -25^L m /5^g IJ +^g u ^L m , is the energy- momentum tensor of bulk matter 
fields, while T nu is the energy-momentum tensor localized on the brane and which is defined 

by T au = — 25 Lbrane/ 8 g 11 " + gfiv Lbrane- 

The delta function 5 (B) denotes the localization of brane contribution. In the 5D space- 
time a brane is a fixed point of the Z2 symmetry. The basic equations on the brane are 
obtained by projections onto the brane world. The induced 4D metric is gjj = ^gu — njnj, 
where nj is the space-like unit vector field normal to the brane hypersurface ^>M. The unit 
vector field satisfies the condition rijn 1 = — 1. In the following we assume ^L m = 0. In the 
brane world models only gravity can probe the extra dimensions. 

Assuming a metric of the form ds 2 = (ninj + gu)dx l dx J , with riidx 1 = d\ the unit 
normal to the \ = constant hypersurfaces and gu the induced metric on \ = constant 
hypersurfaces, the effective 4D gravitational equation on the brane takes the form J4]: 



G 



~Agn V + k\T au + k^Sn V — E 



(3) 



where S nu is the local quadratic energy-momentum correction 



S, 



^TT au - -T«T va + 1^ (32*%, - T' 



(4) 



and E uu is the non-local effect from the free bulk gravitational field, the transmitted pro- 
jection of the bulk Weyl tensor Cjajb, Ejj = CjAjB nAnB , with the property Ejj — > 
Enu^jSj as x ~ ¥ 0. We have also denoted k\ = 8ttG, with G the usual 4D gravitational 
constant. 
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The 4D cosmological constant, A, and the 4D coupling constant, k&, are related by 
A = A;|(A 5 + fcgAf/6)/2 and fcf = respectively. In the limit A^ 1 — > we recover 



a. 



standard general relativity 

The Einstein equation in the bulk and the Codazzi equation also imply the conservation 
of the energy- momentum tensor of the matter on the brane, DJT^ = 0, where D u denotes 
the brane covariant derivative. Moreover, from the contracted Bianchi identities on the 
brane it follows that the projected Weyl tensor obeys the constraint D V E^ = k\D v S^ '. 

The symmetry properties of imply that in general we can decom- 

pose it irreducibly with respect to a chosen 4-velocity field as E^ v = 

-k A U [uyUy + + P^ + 2Q {pL u v) , where k = k 5 /k 4 , = g^ + u^u u projects or- 

thogonal to the "dark radiation" term U = —k~^E tlv u il u v is a scalar, = k~ A h°^E a pvP 
is a spatial vector and P^ u = —k~ 4 a h v ) 13 — ^h^h af3 E a p is a spatial, symmetric and 
trace- free tensor [3|. 

In the case of the vacuum state we have p = p = 0, T^ u = 0, and consequently = 0. 
Therefore the field equation describing a static brane takes the form 

Ryv = —E^v + Ag^, (5) 

with the trace R of the Ricci tensor R^ v satisfying the condition R = R^ = 4A. 

In the vacuum case E^ v satisfies the constraint D V E^ U = 0. In an inertial frame at any 
point on the brane we have u M = Sq and h, w = diag(0, 1, 1, 1). In a static vacuum = 
and the constraint for E^ takes the form 

l -D,U + ^UA, + D V P„ V + A V P^ = 0, (6) 

where = u u D v u^ is the 4-acceleration. In the static spherically symmetric case we 
may chose A^ = A{r)r^ and P^ v = P(r) [r^r v — \h^, where A(r) and P(r) (the "c.ark 
pressure" ) are some scalar functions of the radial distance r, and is a unit radial vector 



12| • 



B. Static and spherically symmetric brane world black holes 

Static black holes are described by the static and spherically symmetric metric given by 

ds 2 = - e ^) dt 2 + e A(r) rfr 2 + r 2 ^2 + sin 2 (?) 
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With the metric given by ([7]) the gravitational field equations and the effective energy 



momentum tensor conservation equation in the vacuum take the form 



1 A'\ 1 

+ ^ = 3aP + A, (8) 



e- x (- + V) - 4 = ol (U + 2P) - A, (9) 



e ->U + ^ + ^-^-)= a{ u-P)-A, (10 ) 



2 V 2 r 2 

U' + 2P' 6P 



v 



(11) 



2P + P r(2P + P)' 
where ' = d/dr, and we have denoted a = WiiG/k^Xi,. 

The field equations (l8l- ffTUl) can be interpreted as describing an anisotropic "matter 
distribution", with the effective energy density p , radial pressure P eS and orthogonal 
pressure Pj^ ff , respectively, so that p cS = 3aU + A, P cfr = all + 2aP — A and Pf 1 = 
all — aP — A, respectively, which gives the condition p eS — P eff — 2P£ ff = 4A = constant. 
This is expected for the 'radiation' like source, for which the projection of the bulk Weyl 
tensor is trace-less, E?. = 0. 

Eq. dSJ) can immediately be integrated to give 

r r 3 

where C\ is an arbitrary constant of integration, and we denoted GM\j (r) = 3a Jq U(r)r 2 dr. 

The function Mu is the gravitational mass corresponding to the dark radiation term 
(the dark mass). For U = the metric coefficient given by Eq. (fT2l must tend to the 
standard general relativistic Schwarzschild metric coefficient, which gives C\ = 2GM, where 
M = constant is the baryonic (usual) mass of the gravitating system. 

By substituting v' given by Eq. (TTTj) into Eq. (jUJ), and with the use of Eq. (fl2"l) . we obtain 
the following system of differential equations satisfied by the dark radiation term U, the 
dark pressure P and the dark mass Mu, describing the vacuum gravitational field, exterior 
to a massive body, in the brane world model 8]: 

dMu ^U. (13) 



dr G 

dU (2U + P) \2GM + GM V - f Ar 3 + a (U + 2P) r 3 



dP 6P 
_ 2 ^--_ , (14) 



dr ~ r 2 (l- - Mll _ A r 2^ dr r 

\ r r 3 



To close the system a supplementary functional relation between one of the unknowns U, 
P and M\j is needed. Generally, this equation of state is given in the form P = P(U). Once 
this relation is known, Eqs. (IT5|) - (jl4p give a full description of the geometrical properties of 
the vacuum on the brane. 

In the following we will restrict our analysis to the case A = 0. Then the system of 
equations (jl~3l) and (fill) can be transformed to an autonomous system of differential equations 
by means of the transformations q = IGMjr + GMjj/r, \i = 3ar 2 U, p = 3ar 2 P, 9 = lnr 
where \i and p are the "reduced" dark radiation and pressure, respectively. With the use of 
the new variables, Eqs. ( |T3l) and (fll]) become 

dq sik\ 

Te=»- q ' (15) 
dfi (2a* + p) [q + W + 2p)] dp 

2 ^ + 2 V-2p- (16) 



d6 1-q d9 

Eqs. (fT3|) and (il4l . or, equivalently, (jl~5l) and (|T6|) . are called the structure equations of 

the vacuum on the brane [8|. In order to close the system of equations f|T5|) and f|T6|) an 

"equation of state" p = relating the reduced dark radiation and the dark pressure 

terms, is needed. Once the equation of state is known, exact vacuum solutions of the 

gravitational field equations on the brane can be obtained. The opposite procedure can 

also be followed, that is, by specifying the functional form of the metric tensor, the dark 

radiation and the dark pressure can be obtained from the field equations. Therefore, several 

exact solutions of the gravitational field equations on the brane can be obtained [12I, Q, Q] • 

C. Rotating brane world black holes 

In order to study rotating black hole solutions on the brane it is convenient to assume 
that the axisymmetric and stationary metric is of the Kerr-Schild form, and can be expressed 
in the form of its linear approximation around the flat metric, ds 2 = (ds 2 ) j lat + H (l^dx^) 2 , 
where is a null, geodesic vector field in both the flat and full metrics, and H is an 
arbitrary scalar function. By introducing a set of coordinates = {u, r, 6, <f)}, the metric 



can be written in the alternative form 



19J 



ds 2 



(du + drf + dr 2 + Y.dd 2 + [r 2 + a 2 ) sin 2 9d(p 2 + 2a sin 2 9drd(f> + H (du — a sin 

(17) 
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where H = H (r,8), the parameter a is related to the angular momentum of the black hole, 
and the quantity E is defined as E = r 2 + a 2 cos 2 9. The condition R = gives for H the 
differential equation 

(MM)— 

with the general solution H = (2Mr — (3) /E, where M and /3 are arbitrary integration 
constants. By applying the Boyer-Lindquist transformation du = dt — (r 2 + a 2 ) dr/A, d(f> = 
dip — adr/A, where A = r 2 + a 2 — 2Mr + f3, the induced metric for a rotating black hole 
on the brane takes the form of the Kerr-Newman solution of standard general relativity, 
describing a stationary and axisymmetric charged black hole. The parameter M can be 
interpreted as the mass of the black hole, but since there is no electric charge on the brane, 
the parameter (3, the tidal charge parameter, which can take both positive and negative 
values, carries the imprints of a non-local, Coulomb type interaction from the bulk. The 
components of the projections of the Weyl tensor from the bulk are given by E\ = —E"£ = 
—(3 [E — 2 (r 2 + a 2 )] /E 3 , E r r = -E e e = /3/E 2 , and = —2/3a (r 2 + a 2 ) sin 2 6>/E 3 , which 
shows a clear analogy with the energy-momentum tensor of a charged rotating black hole in 
standard general relativity jl9 |. 

III. THIN ACCRETION DISKS ONTO BLACK HOLES 

The theory of mass accretion around rotating black holes was developed for the general 



models introduced by Shakura and Sunyaev 



relativistic case by Novikov and Thorne 25 L They extended the steady-state thin disk 



24j to the curved space-time, by adopting the 



equatorial approximation for the stationary and axisymmetric geometry. The time- and 
space-like Killing vector fields (d/dt) >1 and (<9/<90) M describe the symmetry properties of 
this type of space-time, where t and r are the Boyer-Lyndquist time and radial coordinates, 
respectively. 

The horizontal size of the thin disk is negligible as compared to its vertical extension, i.e, 
the disk height H, defined by the maximum half thickness of the disk, is much smaller than 
any characteristic radii r of the disk, H « r. In the steady-state accretion disk models, 
the mass accretion rate Mq is supposed to be constant in time, and the physical quantities 
of the accreting matter are averaged over a characteristic time scale, e.g. At, and over the 
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azimuthal angle A0 = 2tt, for a total period of the orbits and for the height H . The plasma 
moves in Keplerian orbits around the black hole, with a rotational velocity Q, and the plasma 
particles have a specific energy E, and specific angular momentum L, which depend only 
on the radii of the orbits. The particles are orbiting with the four-velocity in a disk 
having an averaged surface density E. The accreting matter is modeled by an anisotropic 
fluid source, where the density po (the specific heat is neglected), the energy flow vector g M 
and the stress tensor t^ v are measured in the averaged rest-frame. The energy-momentum 
tensor describing this source takes the form 

T*" = p uV + 2u ip q v) + e v , 

where = 0, u^ u = 0. The four- vectors of the energy and of the angular momentum 
flux are defined by 

-E* = T» u {d/dt) u J 11 = T" u (d/d(f>y , (19) 

respectively. The four dimensional conservation laws of the rest mass, of the energy and of 
the angular momentum of the plasma provide the structure equations of the thin disk. By 
integrating the equation of the rest mass conservation, V M (po^ M ) = 0, it follows that the 
time averaged accretion rate M is independent of the disk radius: 

M = -27rr£w r = const , (20) 



where a dot represents the derivative with respect to the time coordinate [26( . The averaged 
rest mass density is defined by 

£(r) = I" (f>o)dz, (21) 

where (po) is the rest mass density averaged over At and 2tc. The conservation law V^-E^ = 
of the energy can be written in an integral form as 

[MqE - 27rrflWy], r = Air^f—gFE , (22) 

where a comma denotes the derivative with respect to the radial coordinate r. Eq. fl22|) 
shows the balance between the energy transported by the rest mass flow, the dynamical 
stresses in the disk, and the energy radiated away from the surface of the disk, respectively. 
The torque W^ r in Eq. fl22l is given by 

W/ = f H (t/)dz, (23) 
12 



where (£</, r ) is the — r component of the stress tensor, averaged over At and over a 2tt 
angle. The law of the angular momentum conservation, V^J^ = 0, states in its integral 
form the balance of the three forms of the angular momentum transport, 

[M L - 2nrW/], r = An^g~FL . (24) 



261 ] . This means that 



By eliminating W<f, T from Eqs. (1221) and (|24j) . and by applying the universal energy- 
angular momentum relation dE = QdJ for circular geodesic orbits in the form E >r = QL r , 
the flux of the radiant energy over the disk can be expressed in terms of the specific energy, 
angular momentum and the angular velocity of the black hole. Then the flux integral leads 
to the expression of the energy flux F(r), which is given by 

F(r) = ~ A M '^\ 7 ^ f (E - nl)l r dr , (25) 



4ttv— # (E — QL) 2 - , - 
where the no-torque inner boundary conditions were also prescribed 
the torque vanishes at the inner edge of the disk, since the matter at the marginally stable 
orbit r ms falls freely into the black hole, and cannot exert considerable torque on the disk. 
The latter assumption is valid as long as strong magnetic fields do not exist in the plunging 
region, where matter falls into the hole. 

Once the geometry of the space-time is known, we can derive the time averaged radial 
distribution of photon emission for accretion disks around black holes, and determine the 
efficiency of conversion of the rest mass into outgoing radiation. After obtaining the radial 
dependence of the angular velocity Q, of the specific energy E and of the specific angular 
momentum L of the particles moving on circular orbits around the black holes, respectively, 
we can compute the flux integral (l2"5l) . 



Let us consider an arbitrary stationary and axially symmetric geometry, 

ds 2 = g tt dt 2 + g t<j) dtd(j) + g„dr 2 + geedO 2 + g^dcp 2 , (26) 

where in the equatorial approximation (16* — tt/2\ <C l)the metric functions g u , gt<f>, g rr , Qee 
and g^ depend only on the radial coordinate r. The geodesic equations take the form 

( dt\ 2 _ Eg H + Lg t4> ^ 



dr 



dr 



i 9t<p - 9ttg<t><t> 

* = _Egu t + Lgu^ (2g) 
9t<t> - 9ttg^ 
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and 



9rr(^) =V(r), (29) 



respectively, where r is the affine parameter, and the potential term V(r) is defined by 

V(r) = -1 + ^ + 2EZ^ + ZV (3Q) 

9t<j> ~ 9tt9<t><t> 

For circular orbits in the equatorial plane the following conditions must hold 

V{r) = 0, V r (r) = 0. (31) 

These conditions give the specific energy E, the specific angular momentum L and the 
angular velocity Q of particles moving on circular orbits around spinning general relativistic 

stars as 

E = - , 9tt + 9t ^ , (32) 

-g tt -2g t( p,- " () i 



\J-gtt - 2g t( fl - g^Vl 2 
_ d(f) _ -9t<t>,r + \J(gt<t>,r) 2 - gu,,*, , 



The marginally stable orbit around the central object are determined by the condition 

V rr (r) = 0. (35) 

Let us represent the effective potential in the form V(r) = —1 + f/g, where / = E 2 g^ + 
2ELg tr j> + L 2 g ( j, < j ) and g = g 2 ± — gttg<t><j>, respectively. Then from the condition V(r) = we 
obtain first —1 + f/g = 0, which implies / = g. From V r (r) = we obtain (f\ r g — fg,r) 1 9 2 = 
0, while V rr (r) = gives g~ l {f >rr — g,rr) = 0, since V(r) = and V r (r) = 0. If g ^ we 
have 

E 2 g^^ rr + ZELgt^^r + L 2 g tt<rr — (g 2 ^ — gttg </></>) ,rr = . (36) 

By inserting Eqs. (13"2]) -(133~ |) into Eq. fl36l) . and solving the resulting equation for r, we 
obtain the marginally stable orbits, once the metric coefficients gu, gt<j> and g^ are explicitly 
given. 

In the case of a static and spherically symmetric geometry, given by Eq. ([7]), the geodesic 
equations for particles orbiting in the equatorial plane take the form 



e 2 » 



'£y=*W£)W)=*Wg)^ (3 7 ) 
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and the effective potential can be written as 

V(r) = e?lU ~ ) . (38) 



L 2N 



The conditions for a stable particle orbit are again V(r) = and V r (r) = 0. From these 
conditions we obtain 



>v r e v 



^ = \/^-, (39) 

E ~ Ve» - r 2 Q 2 ' ^ 

L ~ y/e" - r 2 Q 2 ' ^ 
At the marginally stable orbit (or the innermost stable circular orbit) r ms the condition 
V rr (r) = holds, condition from which we can derive the value of r ms for a specified 
function u(r). 

After inserting Eqs. (I3"9l - (j4"ll) into the integral (I25I) . and taking into account that ^fg = 
r 2 e (u+\)/2 £ or q _ n /2, we can compute the flux F(r) over the whole disk surface for each 
brane black hole geometry given by the metric functions z/(r) and A(r). 

The accreting matter in the steady-state thin disk model is supposed to be in thermody- 
namical equilibrium. Therefore the radiation emitted by the disk surface can be considered 
as a perfect black body radiation, where the energy flux is given by F(r) = aT 4 (r) (a is the 
Stefan-Boltzmann constant), and the luminosity L (a;) has a black body spectrum 29]: 



L (u) = Arcd 2 I H = - cos 7 ^ 3 P ^ (42) 

7r J n exp \uj/T) — 1 

Here d is the distance to the source, I(u) is the Planck distribution function, 7 is the 
disk inclination angle, and r» and rj indicate the position of the inner and outer edge of the 
disk, respectively We take = r ms and rj — > 00, since we expect the flux over the disk 
surface vanishes at r — ► 00 for any kind of brane black hole geometry. 

The flux and the emission spectrum of the accretion disks around black holes satisfy 
some simple scaling relations, with respect to the simple scaling transformation of the radial 
coordinate, given by r — * f = r/M, where M is the mass of the black hole. Generally, 
the metric tensor coefficients are invariant with respect of this transformation, while the 
specific energy, the angular momentum and the angular velocity transform as E — > E, 
L — > ML and Q — > Cl/M, respectively. The flux scales as F(r) —>■ F(f)/M A , giving the 
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simple transformation law of the temperature as T(r) — > T (f) /M. By also rescaling the 
frequency of the emitted radiation as u; — > Q = u/M, the luminosity of the disk is given by 
L (u) — > L (u>) /M. On the other hand, the flux is proportional to the accretion rate M , 
and therefore an increase in the accretion rate leads to a linear increase in the radiation 
emission flux from the disc. 

The efficiency e with which the central object converts rest mass into outgoing radiation 
is the other important physical parameter characterizing the properties of the accretion 
disks. The efficiency is defined by the ratio of two rates measured at infinity: the rate of 
the radiation of the energy of the photons escaping from the disk surface to infinity, and 
the rate at which mass-energy is transported to the black hole. If all the emitted photons 
can escape to infinity, the efficiency depends only on the specific energy measured at the 
marginally stable orbit r ms , 

e = 1 - E ms . (43) 

For Schwarzschild black holes the efficiency is about 6%, no matter if we consider the photon 
capture by the black hole, or not. Ignoring the capture of radiation by the black hole, e is 
found to be 42% for rapidly rotating black holes, whereas the efficiency is 40% with photon 
capture in the Kerr potential. 

IV. THIN DISK ACCRETION ONTO BRANE WORLD BLACK HOLES 

In the present Section we consider the accretion properties of several classes of brane world 
black holes, which have been obtained by solving the vacuum gravitational field equations 
for the metrics given by Eqs. (jjT} and ffTT|) . respectively, where the metric functions v and 
A depend only on r. There are many black hole type solutions on the brane, and in the 
following we analyze only four particular example, including three static black hole solutions, 
described by various metric potentials u(r) and A(r), and the rotating generalization of 
the Kerr black hole. All the corresponding metric functions satisfy the gravitational field 
equations on the brane. The energy flux F(r) and the disk emission spectrum is obtained 
for supermassive brane world black holes, with a total mass of 2.5 x 1O 6 M and by using a 
mass accretion rate of 2 x 10 -6 M Q yr _1 . The inclination angle 7 used for the calculation of 
the spectra is set to cos 7 = 1. 
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A. The DMPR brane black hole 



The first brane black hole we consider is a solution of the vacuum field equations, obtained 



by Dadhich, Maartens, Papadopoulos and Rezania in [12|, which represent the simplest 
generalization of the Schwarzschild solution of general relativity. We call this type of brane 
black hole as the DMPR black hole. For this solution the metric tensor components are 
given by 



„ 2M Q 

e v = e~ A = 1 - — + (44) 



where Q is the so-called tidal charge parameter. In the limit Q — > we recover the 
usual general relativistic case. The metric is asymptotically flat, with lim r _ KX) exp (y) = 
lim^ooexp (A) = 1. There are two horizons, given by 



rf = M ± \jM 2 - Q. (45) 

Both horizons lie inside the Schwarzschild horizon r s = 2M, < < r\ < r s . In the 
brane world models there is also the possibility of a negative Q < 0, which leads to only one 
horizon rh+ lying outside the Schwarzschild horizon, 



r h+ = M+ y M 2 + Q > r s . (46) 

In this case the horizon has a greater area than its general relativistic counterpart, so that 
bulk effects act to increase the entropy and decrease the temperature, and to strengthen the 
gravitational field outside the black hole. 

DMPR brane black holes are characterized by the metric functions ( 1441) . If we insert e u 
from these equations into Eq. (1381) . we obtain the effective potential V(r) for this type of 
black hole for any particle with a specific angular momentum L, orbiting around the black 
hole, as a function of the total mass M, and of the tidal charge Q of the black hole. In Fig. [1] 
we present the radial profile of V(r) for M/L = 4, and different values of the tidal charge 
Q, running between — M 2 and M 2 . For comparison we have also plotted the Schwarzschild 
potential, corresponding to Q = 0. 

By increasing Q from zero to M 2 we also increase the potential barrier, as compared to 
the Schwarzschild case, whereas negative tidal charges lowers the barrier, as expected for the 
potential of the Reissner-Nordstrom type black holes. The variation of Q also modifies the 
position of the marginally stable orbit, as shown by the shift of the cut-off in the left hand 
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FIG. 1: The effective potential of a DMPR black hole with total mass M and tidal charge Q for the 
specific angular momentum L = AM. The solid line is the effective potential for a Schwarzschild 
black hole with the same total mass. The values of Q are Q = —M 2 (long dashed line), Q = — 0.5M 2 
(short dashed line), Q = 0.5M 2 (dotted line) and Q = M 2 (dot-dashed line), respectively. 

side of the radial profiles of the photon flux F(r) over the disk surface, which is presented 
in the left plot in Fig. [2j 

For positive increasing tidal charges, r ms takes lower values, and the negative decreasing 
charges lead to the increase of r ms . The position of the maximum of the radiation flux is 
changed similarly due to the variations of r ms . By increasing the value of Q from — M 2 
to M 2 , we shift the radius of the maximal flux to lower and lower values, approaching the 
marginally stable orbit. The vertical shift in the effective potential modifies the amplitude 
of the flux profile as well: higher values of the potential for Q = 0.5M 2 and Q = M 2 
involve a higher specific energy of the orbiting particles, and therefore higher flux values. 
For Q = — 0.5M 2 and Q = — M 2 , the bound state of the particles has a lower energy, and a 
smaller amount of energy is radiated away from the disk surface, which decreases the peak 
in the flux profile. In the right plot of Fig. [2] we present the variation of the flux for different 
values of the mass accretion rate. 

These effects also cause a shift of the cut-off frequency of the black body emission spec- 
trum of the accretion disk, as we can see in the left plot of Fig. |3j The increasing positive 
tidal charges shift the cut-off to higher frequencies, which produces a harder spectrum, 
whereas the negative decreasing values of Q soften the disk spectrum by lowering the cut-off 
frequencies. 
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FIG. 2: Figure on the left: The time-averaged flux radiated by the accretion disk around a DMPR 
black hole for different values of the tidal charge parameter Q: Q = —M 2 (long dashed line), 
Q = -0.5M 2 (short dashed line), Q = 0.5M 2 (dotted line) and Q = M 2 (dot-dashed line), 
respectively. The flux for a Schwarzschild black hole is plotted with a solid line. The mass 
accretion rate is 2 x 10~ 6 MQyr^ 1 . Figure on the right: The time-averaged flux radiated by the 
accretion disk around a DMPR black hole and tidal charge Q = —M 2 for different values of the 
mass accretion rate Mq: 1.0 x lO~ 6 M0yr~ 1 (long dashed line), 1.5 x lO _6 M0yr _1 (short dashed 
line), = 2.0 x lO^Moyr" 1 (solid line), 2.5 x lO^Mey^ 1 (dotted line), and 3.0 x lO^Moyr" 1 
(dot-dashed line), respectively. The total mass of the black hole is 2.5 x 1O 6 M 

The variation of the emission spectrum for different values of the mass accretion rate 
is presented in the right plot in Fig. [3j We also present the conversion efficiency e of the 
mass accreted by the DMPR black hole into radiation measured at infinity, which is given 
by Eq. (|43j) . The value of e measures the efficiency of energy generating mechanism by mass 
accretion. The amount of energy released by the matter leaving the marginally stable orbit, 
and falling into the black hole, is the binding energy E ms of the black hole potential. For 
different tidal charges the values of E ms , together with the radii of the marginally stable 
orbits, are given in Table [J where the figures corresponding to the Schwarzschild black hole 
corresponds to Q = 0. As we increase Q from — M 2 to M 2 , r ms is decreasing from values 
greater than the radius of the marginally stable orbit for the Schwarzschild geometry to 
lower ones. The efficiency has an opposite trend: for negative tidal charges it has smaller 
values than in the case of the Schwarzschild black holes, and it exceeds the latter for Q > 0. 
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FIG. 3: Figure on the left: The emission spectrum of the accretion disk around a DMPR black hole 
for different values of the tidal charge parameter Q: Q = —M 2 (long dashed line), Q = — 0.5M 2 
(short dashed line), Q = 0.5M 2 (dotted line) and Q = M 2 (dot-dashed line), respectively. The 
solid line represents the disk spectrum for a Schwarzschild black hole. The mass accretion rate 
is 2 x 10 _6 M Q yr~ 1 . Figure on the right: The emission spectrum of the accretion disk around a 
DMPR black hole with tidal charge Q = —M 2 for different values of the mass accretion rate Mq: 
1.0 x lO~ 6 M yr~ 1 (long dashed line), 1.5 x 10~ 6 M o yr- 1 (short dashed line), 2.0 x 10~ 6 M o yr~ 1 
(solid line), 2.5 x 10~ 6 M Q yr _1 (dotted line), and 3.0 x lO _6 M yr _1 (dot-dashed line), respectively. 
The total mass of the black hole is 2.5 x 1O 6 M 



Q [M 2 ] 


Tms [M] 


e 


-1 


7.3100 


0.0476 


-0.5 


6.6949 


0.0517 





6.0000 


0.0572 


0.5 


5.1695 


0.0655 


1 


4.0019 


0.0814 



TABLE I: The marginally stable orbit and the efficiency for different DMPR black hole geometries. 
The case Q = corresponds to the standard general relativistic Schwarzschild black hole. 

B. The CFM brane black hole 



Two families of analytic solutions in the brane world model, parameterized by the ADM 
mass and the PPN parameters (3 and 7, and which reduce to the Schwarzschild black hole 
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for (3 = 1, have been found by Casadio, Fabbri and Mazzacurati in [14j. We call the 
corresponding brane black holes as the CFM black holes. 
The first class of solutions is given by 

2M 



1 - 



and 



i 3M 



4 

9^ 



(47) 



(48) 



1 - 

respectively, where 77 = 7 — l = 2(/3 — 1). As in the Schwarzschild case the event horizon 
is located at r = = 2M. The solution is asymptotically flat, that is lim^ooe^ = e Uao = 
lim^oo e A = e A °° = 1 . 

The second class of solutions corresponding to brane world black holes obtained in [ijj] 
has the metric tensor components given by 

2 

(49) 



77 + ^/1-^(1 + 77) 
1 + 77 



and 



2M 



I + 77) 



(50) 



respectively. The metric is asymptotically flat. In the case 77 > 0, the only singularity in the 
metric is at r = r$ = 2M (1 + 77), where all the curvature invariants are regular, r = ro is a 
turning point for all physical curves. For 77 < the metric is singular at r = = 2M/ (1 — 77) 
and at r , with > r . defines the event horizon. 

The first and second class of the CFM brane black holes are characterized by the metric 
potentials ( 1471) and ( )49l) . Since for the first class of solutions the metric function e u coincides 
with the one of the Schwarzschild black hole, their effective potentials V(r) are the same 
for equal total masses and fixed L. As a consequence, the specific energy, specific angular 
momentum and angular velocity of the particles orbiting around the first class of CFM 
black holes are equal to those of the particles moving at the same Keplerian orbit in the 
Schwarzschild potential. 

However, the radiation flux from the accretion disk shows considerable differences as 
compared to the Schwarzschild case when the parameter 77 is varied. This behavior, shown 
in the left plot of Fig. HI is due to the fact that the proper volume used in the calculation of 
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any integral in this spacetime depends on both the metric functions v{r) and A(r), and in 
the flux integral given by Eq. ([25]) we have 



Since the left hand side of Eq. fl5Tl) is decreasing for negative values of 77, we obtain higher 
and higher flux values by decreasing this metric parameter for 77 < 0. We have the opposite 
effect for 77 > 0, where the increase of this parameter causes the proper volume to also 
increase, and, in turn, the amplitudes of the flux profile to decrease. Because the effective 
potential does not depend on 77, there is no variation in the radius of the marginally stable 
orbit. However, a cut-off appears in the left hand side of the flux profiles for 77 > 0, where 
the denominator in Eq. (I5TI) becomes negative. This gives the criterion 



As one can see in left plot of Fig. [4j for 77 = 5 this criterion is true for any radii greater than 
r ms , but we obtain a cut-off in the left hand side of the flux profile by setting 77 equal to 
77 = 10. 

Although the position of the marginally stable orbit does not change for different values 
of 77, the maximum of the radial distribution of the flux is located at higher and higher radii 
as we increase 77. The right plot in Fig. H] presents the variation of the flux for a fixed mass 
and parameter 7/, and for different values of the accretion rate. 

The variation in the shape and the amplitude of the flux profile for different values of 77 
has a clear effect of shifting the cut-off frequency in the disk spectra. The left plot of Fig. [5] 
shows that the cut-off value shifts to higher frequencies when rj is negative, hardening the 
spectra. For 77 > the disk spectrum becomes softer, with lower cut-off frequencies. The 
effect of the variation of the mass accretion rate on the emission spectra for a fixed mass 
and rj is presented in the right plot of Fig. [51 

Since the effective potential of the CFM black holes of the first class is the same as the 
one of the Schwarzschild geometry, their conversion efficiencies e of the mass accreted by the 
black hole into radiation measured at infinity are also the same. 

If we consider the second class of the CFM black holes, given by Eqs. ( pf9l) and fl50l) . 
respectively, the variation of the parameter 77 causes similar effects in the behavior of the 
photon flux emitted by the accretion disk and its spectrum. However, the metric function 




.1/2 



(51) 




(52) 



22 





FIG. 4: Figure on the left: The time-averaged flux radiated by an accretion disk around a first 
class CFM black holes for various values of rj: 77 = — 10 (long dashed line), rj = — 5 (short dashed 
line), rj = 5 (dotted line) and 77 = 10 (dot-dashed line). The flux for a Schwarzschild black hole 
is plotted with a solid line. The mass accretion rate is 2 x 1O~ 6 M yr" 1 . Figure on the right: 
The time-averaged flux radiated by the accretion disk around a first class CFM black hole for 
■q = —10 and for different values of the mass accretion rate Mq: 1.0 x lO _6 M0yr _1 (long dashed 
line), 1.5 x lO" 6 Af yr" 1 (short dashed line), 2.0 x lO^Moyr" 1 (solid line), 2.5 x lO" 6 M yr" 1 
(dotted line), and 3.0 x lO _6 M0yr~ 1 (dot-dashed line), respectively. The total mass of the black 
hole is 2.5 x 1O 6 M . 

u(r) now differs from the Schwarzschild black hole case, and we obtain different effective 
potentials for different values of the 77. The radial profiles of V(r) are shown in Fig. 0, where 
77 is set to values between —0.8 and 0.8. 

By decreasing the values of 77 from zero to small negative values, we can increase the 
potential barrier around the black hole, and decrease the radius of the marginally stable 
orbit. The minimum of the effective potential is also increased. Any increment in 77 from 
zero to small positive values causes the opposite effects for V(r). Then the parameters 77 < 
increase the flux values and shift r ms to lower radii, while negative values of 77 give lower 
fluxes with a cut-off at higher values of r ms , as shown in the left plot of Fig. [7J 

The variation of the numerical value of the radial coordinate r at the position where 
the radiation flux takes its maximum value follows the tendency of r ms : with increasing rj 
we obtain higher and higher orbits for F max . These effects result in similar shifts in the 
disk spectrum as in the case of the CFM black holes of the first class: in the left plot of 
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FIG. 5: Figure on the left: The emission spectrum of the accretion disk around a first class CFM 
black hole for different values of rj: r/ = — 10 (long dashed line), r] = —5 (short dashed line), 
T) = 5 (dotted line) and 77 = 10 (dot-dashed line). The spectrum for a Schwarzschild black hole 
is plotted with a solid line. The mass accretion rate is 2 x 1O~ 6 M0 yr" 1 . Figure on the right: 
The emission spectrum of the accretion disk around a first class CFM black hole for 77 = — 10 
and for different values of the mass accretion rate Mq\ = 1.0 x lO~ 6 M0yr -1 (long dashed line), 
1.5 x lO" 6 M yr- 1 (short dashed line), = 2.0 x 10' 6 M o yr~ 1 (solid line), 2.5 x lO^Mgyr -1 (dotted 
line), and 3.0 x lO _6 M0yr _1 (dot-dashed line), respectively. The total mass of the black hole is 
2.5 x 1O 6 M . 

Fig. [H] negative values of 77 harden the spectrum, by shifting the frequency cut-off to higher 
frequencies, whereas the spectra become softer, with lower cut-off values, for increasing 
positive values of the metric parameter. 

The right plots of Figs. [7] and [8] show the variation of the flux and emission spectra for a 
fixed mass and 77 and for different values of the accretion rate. 

We present the marginally stable orbits r ms and the conversion efficiency e of the matter 
accreted by the CFM black holes of the second class into radiation measured at infinity in 
Table HT1 With the increasing values of the parameter rj, the radius of the marginally stable 
orbit is also increasing, as already seen in the flux diagram, while e drops down from values 
greater than 6 to lower ones. 
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FIG. 6: The effective potential of the thin accretion disk around a CFM black hole of the second 
class with a total mass of M and L = AM for different values of the parameter 77. The potential 
V(r) corresponding to the Schwarzschild black hole is plotted with a solid line. The values of 77 
are 77 = —0.8 (long dashed line), 77 = —0.4 (short dashed line), 77 = 0.4 (dotted line) and 77 = 0.8 
(dot-dashed line), respectively. 



V 


r ms [M] 


e 


-0.8 


4.3445 


0.0795 


-0.4 


5.1165 


0.0657 





6.0000 


0.0572 


0.4 


6.9722 


0.0501 


0.8 


8.0091 


0.0444 



TABLE II: The marginally stable orbit and the efficiency for CFM black holes of the second kind. 
The case 77 = corresponds to the standard general relativistic Schwarzschild black hole. 

C. The BMD brane world black hole 



Several classes of brane world black hole solutions have been obtained by Bronnikov, 
Melnikov and Dehnen in 18j (for short the BMD black holes). In the following we analyze 



the accretion properties of a particular class of these models, with metric given by 

I - — I . < A = ( I - — I . (53) 
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FIG. 7: Figure on the left: The time-averaged flux radiated by an accretion disk around a second 
class CFM black hole for different values of r/: rj = —0.8 (long dashed line), r] = —0.4 (short dashed 
line), rj = 0.4 (dotted line) and r] = 0.8 (dot-dashed line). The flux for a Schwarzschild black hole 
is plotted with a solid line. The mass accretion rate is 2 x 1O _6 M yr^ 1 . Figure on the right: The 
time-averaged flux radiated by the accretion disk around a second class CFM black hole with the 
parameter rj = 0.8 for different values of the mass accretion rate Mq\ l.Ox lO _6 M yr _1 (long dashed 
line), 1.5 x lO^Moyr -1 (short dashed line), 2.0 x lO" 6 M yr~ 1 (solid line), 2.5 x lO^Msyr -1 
(dotted line), and 3.0 x lO _6 M yr _1 (dot-dashed line), respectively. The total mass of the black 
hole is 2.5 x 1O 6 M . 

where s GN. The metric is asymptotically flat, and at r = — 2M these solutions have a 
double horizon. 

Eqs. (153|) determine the geometry of the BMD brane black holes, which have the effective 
potential plotted in Fig. [9j In the figure we have plotted V(r) for values of s between s = 5 
and s = 20. 

With increasing s the potential barrier increases as well, but we obtain smaller and smaller 
radii for the marginally stable orbits. Although the higher values of s increase the potential 
over the region of the stable Keplerian orbits, and the energy of the orbiting particles, as 
compared to the case of the Schwarzschild potential, the value of # used to calculate the 
flux integral increases more rapidly. Therefore Eq. (|25|) gives smaller fluxes for higher values 
of s. This effect is shown in the left plot of Fig. [TOj where we present the plots of the photon 
flux emitted by the accretion disk in the BMD brane black hole geometry. 

The relative shift of r ms to lower orbits for increasing s can also be well studied in the 
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FIG. 8: Figure on the left: The emission spectrum of the thin accretion disk around a second class 
CFM black hole for different values of r\: r\ = —0.8 (long dashed line), rj = —0.4 (short dashed 
line), i] = 0.4 (dotted line), and rj = 0.8 (dot-dashed line). The spectrum for a Schwarzschild black 
hole is plotted with a solid line. The mass accretion rate is 2 x 1O~ 6 M yr _1 . Figure on the right: 
The emission spectrum of the accretion disk around a second class CFM black hole with parameter 
rj = 10 for different values of the mass accretion rate M$: 1.0 x 10 _6 M Q yr _1 (long dashed line), 
1.5 x 10- 6 M Q yr~ 1 (short dashed line), 2.0 x lO- 6 M yr" 1 (solid line), 2.5 x 10- 6 M Q yr~ 1 (dotted 
line), and 3.0 x lO _6 M yr _1 (dot-dashed line), respectively. The total mass of the black hole is 
2.5 x 1O 6 M . 
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FIG. 9: The effective potential for BMD black holes of a fixed total mass M for L = 4M and 
different values of s. The solid line is the effective potential for a Schwarzschild black hole with 
the same total mass M. The parameter s is set to s = 5 (long dashed line), s = 10 (short dashed 
line), s = 15 (dotted line) and s = 20 (dot-dashed line). 
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FIG. 10: Figure on the left: The time-averaged flux radiated by an accretion disk around a BMD 
black hole for different values of s: s = 5 (long dashed line), s = 10 (short dashed line), s = 15 
(dotted line) and s = 20, respectively. The flux for a Schwarzschild black hole is plotted with a 
solid line. The mass accretion rate is 2 x 1O _6 M0 yr -1 . Figure on the right: The time-averaged flux 
radiated by the accretion disk around a BMD black hole with s = 20 for different values of the mass 
accretion rate Mq: 1.0 x 10~ 6 MQyr _1 (long dashed line), 1.5 x lO _6 M0yr _1 (short dashed line), 
2.0 x lO^Moyr" 1 (solid line), 2.5 x lO^M^r" 1 (dotted line), and 3.0 x lO _6 M yr -1 (dot-dashed 
line), respectively. The total mass of the black hole is 2.5 x 10 6 Mq. 

plot. The maximum of the flux value has the same shift: by increasing s the maximum of 
F(r) is obtained at lower and lower radii. This behavior of the radiation flux results in the 
softening of the emission spectrum of the disk. As seen in the left plot of Fig. [TTJ the cut-off 
of the disk spectrum shifts to lower frequencies for increasing values of s, as compared to 
the ones obtained in the case of the Schwarzschild geometry. The right plots in Figs. [10] and 
[TTl show the effect of the change of the accretion rate for a BMD brane world black hole for 
a fixed mass an s. 

The marginally stable orbit and the conversion efficiency e of the accreted mass into 
radiation measured at infinity for the BMD black holes are presented in Table IIHI Both 
r ms and e have values less than those derived in the Schwarzschild geometry, and they exhibit 
a decreasing tendency as we increase the values of the parameter s. 
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FIG. 11: Figure on the left: The emission spectrum of the thin accretion disk around a BMD black 
hole for different values of the parameter s: s = 5 (long dashed line), s = 10 (short dashed line), 
s = 15 (dotted line), and s = 20. The solid line represents the spectrum for the Schwarzschild black 
hole. The mass accretion rate is 2 x 10~ 6 Mq yr" 1 . Figure on the right: The emission spectrum 
of the accretion disk around a BMD black hole with parameter s = 20 for different values of the 
mass accretion rate Mq: 1.0 x lO _6 M0yr _1 (long dashed line), 1.5 x 10~ 6 MQyr _1 (short dashed 
line), = 2.0 x 10~ 6 M o yr^ 1 (solid line), 2.5 x 10" 6 M Q yr" 1 (dotted line), and 3.0 x 10~ 6 M o yr" 1 
(dot-dashed line), respectively. The total mass of the black hole is 2.5 x 10 6 Mq. 



s 


r ms [M] 


e 


5 


4.8005 


0.0303 


10 


4.4007 


0.0170 


15 


4.2708 


0.0118 


20 


4.2046 


0.0090 



TABLE III: The marginally stable orbit and the efficiency for BMD black holes. 
D. The AG rotating brane world black hole 



The generalization of the Kerr solution of standard general relativity for rotating brane 
world black holes was obtained by Aliev and Gumrukcuoglu [19J, and is given by 



ds 2 



2Mr - /3 N 



2a(2Mr-(3) . 



I - ! r// 2 - -" siir Mills + ^dr 2 ^<!0 2 
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(r* + a 2 + ^ V sin 2 j sin 2 (54) 

We call this solution the AG rotating brane world black hole. The event horizon of the 
black hole is determined by the solution of the equation A = 0, with the largest root given 
by r + = M + \JM 2 — a 2 — (3. The event horizon does exist if the condition M 2 > a 2 + (3 is 
fulfilled. For a negative tidal charge, as a — > M, r+ — > M + y/—(3 > M, a condition that 
is not allowed in standard general relativity. On the other hand for a negative tidal charge 
the extreme horizon r + = M corresponds to a black hole with rotation parameter a greater 
than its mass M. 

Near to the equatorial plane, \9 — 7r/2| <C 1, by introducing the coordinate z = rcosO ~ 
r{9 — 7r/2), the approximate form for the geometry of the AG rotating brane black hole can 
be written as 

ds 2 = -9^~ l dt 2 + r 2 £/{d<j) - udt) 2 + 9~ l dr 2 + dz 2 (55) 
with the metric functions 

si = 1 + a 2 Jrl + 2a 2 Jrl - a%/rj , (56) 
0=1- 2/r* + a 2 Jrl - ft/r^ (57) 
w = 2Mar -3 ^ -1 , (58) 

where we denoted a* = a/M, = r/M, and /3* = /3/M 2 , respectively. The effective 
potential per unit mass for the radial motion is given by 

V(r) = [(r 2 + a 2 ) r 2 + (2Mr - (3) a 2 ] E 2 - (r 2 - 2Mr + pi) L 2 - 2a (2Mr — (3) EL — r 2 A. 

(59) 

The variation of the potential V(r) is represented, as a function of r/M, and for different 
values of the tidal charge parameter /3, in Fig. [121 By increasing f3 from zero to 2M 2 we 
also increase the potential barrier as compared with the standard general relativistic Kerr 
black hole case, whereas negative tidal charges lowers the potential barrier. The changes in 
the value of (3 also modify the positions of the marginally stable orbits. 

In Fig. [13] we present the flux profiles of the accretion disk in the modified Kerr geometry 
( I55p as a function of r/M and for different values of the tidal charge j3 (left plot) and of 
the accretion rate M (right plot). The variation of the numerical value of the tidal charge 
determines similar modifications for the flux values as in the case of the accretion disk around 
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FIG. 12: The effective potential of a rotating AG brane world black hole of a total mass 2.5 x 1O 6 M0 
and with spin a = 0.9982 for the specific energy E = 0.8 and the specific angular momentum 
L = AM. The solid line is the effective potential for a rotating Kerr black hole (/? = 0) with the 
same total mass. The different values of j3 are f3 = —2M 2 (long dashed line), = —M 2 (short 
dashed line), (5 = M 2 (dotted line) and j3 = 2M 2 (dot-dashed line), respectively. 

the DMPR black holes, which can be considered as the static limit (a = 0) of the rotating 
AG black hole. The left hand plots in Figs. [T3]and[2J respectively, show the same variation 
of the flux profiles as a function of the tidal charge. We note that since Q corresponds to 
—(3, an increase in the numerical values of (3 from negative values to positive ones decreases 
the magnitude of the flux, and increases the radius of the marginally stable orbits. The right 
hand plots in Figs. [2] and [13], respectively, also exhibit the same tendency: for higher mass 
accretion rates the flux will be amplified as well. In Figs. [2] and H3] the cut-off values of the 
spectra decrease with the increasing values of (3 and —Q, respectively. The same analogy is 
valid for disk spectra for the static and the rotating cases. The emission spectra in the case 
of the rotating AG brane world black holes are represented, for different values of the spin 
parameter and of the accretion rate, in Figs. [HI 

By comparing Tables [I] and [IV] we can see the same effect of the variation of the tidal 
charge on the efficiency e for the static DMPR and for the rotating AG brane black hole, 
respectively. As (3 and — Q increase from negative values to positive ones, the marginally 
stable orbits shift to higher radii, and the efficiency of the conversion of the accreting mass 
to radiant energy decreases, from values higher than 0.3241 (the efficiency for the standard 
Kerr black hole with a* = 0.9982) to lower ones. 
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FIG. 13: Figure on the left: The time-averaged flux radiated by the accretion disk around a rotating 
AG brane world black hole with spin a = 0.8892 for different values of the tidal charge parameter (3: 
(3 = -2x 10" 3 M 2 (long dashed line), (3 = -1CT 3 M 2 (short dashed line), = 1(T 3 M 2 (dotted line), 
and (3 = 2 x 1CT 3 M 2 (dot-dashed line), respectively. The flux for a Kerr black hole with the same 
total mass and spin is plotted with a solid line. The mass accretion rate is 2 x lO _6 M yr~ 1 . Figure 
on the right: The time- averaged flux radiated by the accretion disk around a rotating AG brane 
world black hole with spin a = 0.8892 and tidal charge (3 = — 2 x 10~ 3 M 2 for different values of the 
mass accretion rate Mq: 1.0 x 10 _6 MQyr _1 (long dashed line), 1.5 x lO _6 M0yr" 1 (short dashed 
line), = 2.0 x lO^M^r" 1 (solid line), 2.5 x lO^M^r" 1 (dotted line), and 3.0 x lO^Moyr" 1 
(dot-dashed line), respectively. The total mass of the black hole is 2.5 x 1O 6 M . 



(3 [10~ 3 M 2 ] 


r ms [M] 


e 


-2 


1.1677 


0.3449 


-1 


1.2019 


0.3329 





1.2277 


0.3241 


1 


1.2511 


0.3169 


2 


1.2716 


0.3109 



TABLE IV: The marginally stable orbit and the efficiency for different rotating AG black hole 
geometries for a* = 0.9982. The case (3 = corresponds to the standard general relativistic Kerr 
black hole. 
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FIG. 14: Figure on the left: The emission spectrum of the accretion disk around a rotating AG 
brane world black hole with spin a = 0.8892 for different values of the tidal charge parameter: 
(3 = -2 x 1(T 3 M 2 (long dashed line), (3 = -1(T 3 M 2 (short dashed line), (3 = 1CT 3 M 2 (dotted 
line), and j3 = 2 x 10~ 3 M 2 (dot-dashed line), respectively. The flux for a Kerr black hole with the 
same total mass and spin is plotted with a solid line. The mass accretion rate is 2 x 10 _6 M Q yr _1 . 
Figure on the right: The emission spectrum of the accretion disk around a rotating AG brane 
world black hole with spin a = 0.8892 and tidal charge (3 = — 2 x 10~ 3 M 2 for different values of the 
mass accretion rate Mq: 1.0 x 10 M0yr _1 (long dashed line), 1.5 x lO _6 M0yr~ 1 (short dashed 
line), = 2.0 x lO^M^r" 1 (solid line), 2.5 x lO^Afoyr -1 (dotted line), and 3.0 x lO^Moyr" 1 
(dot-dashed line), respectively. The total mass of the black hole is 2.5 x 10 6 M Q . 

V. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have considered the basic physical properties of matter forming 
a thin accretion disk in the space-time metric of the brane world black holes. The physical 
parameters of the disk-effective potential, flux and emission spectrum profiles have been 
explicitly obtained for several classes of black holes, and for several values of the parameters 
characterizing the vacuum solution of the generalized field equations in the brane world 
models. All the astrophysical quantities related to the observable properties of the accretion 
disk can be obtained from the black hole metric. 

There are many effective 4D solutions of the vacuum field equations on the brane, with 
arbitrary parameters which depend on properties of the bulk, or are simply put in by using 
general physical considerations. At the present moment it is theoretically not known whether 
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these parameters should be universal over all brane world black holes, or whether each 
separate black hole may have different values of them. Conversely, there is not a single 
complete solution, in the sense that the metric in the bulk is uniquely known. This situation 
is unsatisfactory from a purely theoretical point of view, and a solution of this problem 
seems to be very difficult to be found. Therefore it may be useful to solve the problem of 
the existence and nature of the brane world black holes by investigating more closely the 
existing observational evidence of the black hole properties, and try to discriminate between 
different black hole models by using the data provided by the observational study of the 
astrophysical processes around black holes. 

Testing strong field gravity and the detections of the possible deviations from standard 
general relativity, signaling the presence of new physics, remains one of the most important 
objectives of observational astronomy. Due to their compact nature, black holes provide an 
ideal environment to do this. Presently, the best constraints on the brane world black hole 
parameters can be obtained from the classical tests of general relativity (perihelion preces- 
sion, deflection of light, and the radar echo delay, respectively). The existing observational 
solar system data on the perihelion shift of Mercury, on the light bending around the Sun 
(obtained using long-baseline radio interferometry) , and ranging to Mars using the Viking 
lander, as applied to the DMPR black hole, can constrain the numerical values of both 



the bulk tidal parameter Q and of the brane tension |32j. The stronger limit is obtained 
from the perihelion precession, \Q\ < 6 x 10 7 — 5 x 10 8 cm 2 . An improvement of one order 
of magnitude in the observational data on Mercury's perihelion shift could provide a very 
precise estimate of the bulk tidal parameter, as well as of the brane tension A. 

Observations in the near-infrared (NIR) or X-ray bands have provided important infor- 
mation about the spin of the black holes, or the absence of a surface in stellar type black hole 
candidates. In the case of the source Sgr A*, where the putative thermal emission due to the 
small accretion rate peaks in the near infrared, the results are particularly robust. However, 
up to now, these results have confirmed the predictions of the general relativity mainly in 
a qualitative way, and the observational precision achieved cannot distinguish between the 
different proposed theories of gravitation. However, important technological developments 
may allow to image black holes directly [33|]. A background illuminated black hole will 
appear in a silhouette with radius \^27GM/c 2 , with an angular size of roughly twice that 
of the horizon, and may be directly observed. With an expected resolution of 20 //as, sub- 
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millimeter very-long baseline interferometry (VLBI) would be able to image the silhouette 
cast upon the accretion flow of Sgr A*, with an angular size of ~ 50 //as, or M87, with an 
angular size of ~ 25 //as. For a black hole embedded in an accretion flow, the silhouette will 
generally be asymmetric regardless of the spin of the black hole. Even in an optically thin 
accretion flow asymmetry will result from special relativistic effects (aberration and Doppler 
shifting). In principle, detailed measurements of the size and shape of the silhouette could 
yield information about the mass and spin of the central black hole, and provide invaluable 
information on the nature of the accretion flows in low luminosity galactic nuclei. 

Due to the differences in the space-time structure, the brane world black holes present 
some important differences with respect to their disc accretion properties, as compared to 
the standard general relativistic Schwarzschild and Kerr cases. Therefore, the study of 
the accretion processes by compact objects is a powerful indicator of their physical nature. 
Since the conversion efficiency, as well as the flux and the spectrum of the black body 
radiation in the case of the brane world black holes is different as compared to the standard 
general relativistic case, the astrophysical determination of these physical quantities could 
discriminate, at least in principle, between the different gravity theories, and give some 
constrains on the existence of the extra dimensions. 
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